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We study the time until a filter becomes clogged due to the trapping of suspended particles as 
they pass through a porous medium. This trapping progressively impedes and eventually stops the 
flow of the carrier fluid. We develop a simple description for the pore geometry and the motion of 
the suspended particles which, together with extreme-value statistics, predicts that the distribution 
of times until a filter clogs has a power-law long-time tail, with an infinite mean clogging time. These 
results and its consequences are in accord with simulations on a square lattice porous network. 
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In this letter, we investigate the time required for a 
filter to clog. In a typical filtration process, a dirty fluid 
is "cleaned" by passing it through a porous medium to 
remove the suspended particles. The medium enhances 
filtering efficiency by increasing both the available filter 
surface area for trapping suspended particles, as well as 
the exposure time of the suspension to the active surfaces. 
Such a mechanism is the basis of water purification, air 
filtration, and many other separation processes |l|,|| . As 
suspended particles become trapped, the fluid permeabil- 
ity of the medium gradually decreases, and eventually the 
filter becomes clogged. Determining the time dependence 
of this clogging is basic to predicting when a filter is no 
longer useful, either because of reduced throughput or 
reduced filtration efficiency, and should be discarded. 

We develop a minimalist model which provides an in- 
tuitive understanding for the clogging of a filter. Our 
analytical results are based on representing the clogging 
of a porous medium by the clogging of a single parallel 
array of pores which are blocked in decreasing size order 
(Fig. |l]) . This geometrical reduction is based on previous 
filtration studies which showed that, when particles and 
pores are of comparable size, clogging preferentially oc- 
curs in the upstream end of the network J^B]. A single 
parallel array of pores represents the ultimate limit of 
this gradient-controlled process. 

The approximation of size-ordered blocking is based on 
the commonly-used picture that a particle has a proba- 
bility proportional to the relative flux to enter a given 
unblocked pore from among the outgoing pores at a junc- 
tion [^U^H- For Poiseuille flow, this flux is proportional 
to r 4 , where r is the pore radius. Size-ordered blocking 
arises if the exponent 4 in this flow- induced pore entrance 
probability is replaced by oo. For a broad distribution of 
pore radii, this is not a drastic approximation. This or- 
dering also provides a direct correspondence between the 
radius of the currently-blocked pore and the time until 
it is blocked. We then use extreme value statistics || 
to compute the radius distribution of the last few un- 
blocked pores, together with the connection between the 
radius of the currently blocked pore and the network per- 
meability, to determine the distribution of filter clogging 
times. The predictions of this simple modelling are in 



good agreement with Monte Carlo simulations of clog 
ging on a square lattice porous network. 
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FIG. 1. Illustration of the clogging evolution in a parallel 
array of w bonds of distributed radii. Bonds are blocked in 
decreasing size order. The overall flow decreases more slowly 
in the latter stages (vertical arrow). The total flow is in- 
versely related to the time until each blockage event, tk, with 
k = w, w — 1, . . . ,. 

In our simulations of clogging, we consider a lattice 
porous medium whose bonds represent pores. We assume 
Poiseuille flow, in which the fluid flux through a bond of 
radius r, is proportional to — rf Vp, where Vp is the lo- 
cal pressure gradient when a fixed overall pressure drop 
is imposed. Dynamically neutral suspended particles tra- 
verse the medium in accordance with the local flow and 
the flow-induced entrance probability at each junction. A 
crucial feature is that particles are injected at a finite rate 
which is proportional to the overall fluid flux, correspond- 
ing to a fixed non-zero density of suspended particles in 
the fluid || . This is in distinction to many previous stud- 
ies of filtration, where particles were injected singly and 
tracked until trapping occurred [§[j4); this corresponds 
to an arbitrarily dilute suspension. While useful for un- 
derstanding the percolation process induced by clogging 
|}|-|5],f7|, the time evolution of filtration necessitates the 
consideration of a finite-density suspension. 
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For the trapping mechanism we adopt size exclusion, 
in which a particle of radius r par ti c ie is trapped within 
the first bond encountered with rb on d < ^particle f|- We 
assume that a trapped particle blocks a pore completely 
and permanently. The difference between partial and 
total blockage in a single trapping event appears to be 
immaterial for long-time properties (J] . The overall trap- 
ping rate is then controlled by the relative sizes of parti- 
cles and bonds. For simplicity and because it occurs in 
many types of porous media H, we consider the Hertz 
distribution of pore and particle radii, respectively, 



b(r) = 2are 



and 



p(r) = 2re 



(1) 



with the ratio between the average bond and particle 
radii s = 1/ i/a a basic parameter which determines the 
nature of the clogging process. 

First consider the case where pores are larger than par- 
ticles (s > 1). Many particles must be injected before 
a sufficiently large particle arises which can block the 
largest pore. Since the permeability of the system re- 
mains nearly constant when only a few pores are blocked, 
the time between successive particle injection events dur- 
ing this initial stage is nearly constant. Once the largest 
pores are blocked, it takes many fewer particles to block 
the remaining smaller pores and clogging proceeds more 
quickly. Thus the clogging time is dominated by the ini- 
tial blockage events. 

To estimate the clogging time in this large-pore limit, 
let us find the number of injection events before encoun- 
tering a particle large enough to block the largest bond. 
The radius of this bond is given by the criterion H 
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that one bond out of w has has radius r max or larger. 
This gives r^ix = s-\/ln w. Similarly, the largest parti- 
cle out of N has radius rmax = VhiiV. Thus N± rj w s2 
particles typically need to be injected before one occurs 
which is large enough to block the largest bond. Continu- 
ing this picture sequentially, the radius of the fc th -largest 
bond is given by Eq. (||), but with l/w replaced by k/w; 
this gives r^ = s^J\n{w/k). Therefore N k ps (w/kf 2 
particles need to be injected before the /c th -largest bond 
is blocked. Since Nk decreases rapidly with k, the ini- 
tial blockage events control the clogging time T, whose 
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lower bound is given by T sa Nx/w > w s . The factor 
of l/w arises because the particle injection rate is pro- 
portional to w for a fixed-density suspension and a fixed 
pressure drop. However, there is considerable particle 
penetration in the large-pore limit, so that the equiva- 
lence to the clogging of a parallel bond array is not really 
appropriate and the bound on T is quite crude. 

The latter case where pores are smaller than particles 
(s < 1) is simpler, as each particle injection event typi- 
cally leads to the clogging of the first pore entered. Now 
the initial pores are blocked quickly, while later pores 



are blocked more slowly because the overall flow rate de- 
creases significantly near the end of the clogging process. 
We argue that clogging is dominated by the time of these 
later blockage events. First, let us determine how the 
flow rate varies as the last few bonds get blocked. Since 
only the smallest bonds remain open near clogging, the 
permeability is determined by these smallest radii. We 
estimate the radius of the fc th smallest bond from 



2are- ar dr = -, 



(3) 



which gives r k = s ^Jk/w. The permeability of a parallel 
bundle of these k smallest pores is then 
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k 
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while the initial system permeability (obtained by setting 
k = w above) is simply k(w) = s w. 

Since the overall fluid flow is proportional to the per- 
meability for a fixed pressure drop, the time increment tk 
between blocking the (k — l) st -smallest and fc th -smallest 
pore scales as 



tk 



wn{k) 



w 
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Here the factor w in the denominator again accounts for 
a particle injection rate proportional to w. The clogging 
time T is now dominated by the time to block the small- 
est bond, so that T > t\ oc w 2 . 

Thus the relative pore and particle radii drives a tran- 
sition in the clogging process. For s 2 > 3, clogging is 
dominated by the initial blockage events and T > w s , 
while for s 2 < 3, the last events dominate, leading to 
T > w 2 . For the small-pore limit, we may carry the 
analysis further and obtain the distribution of clogging 
times. This distribution gives a divergent mean clogging 
time; nevertheless, suitably defined moments of the clog- 
ging time distribution scale as w 2 . 

We find the clogging time distribution in terms of the 
radius distribution of the smallest bond, since this bond 
ultimately controls clogging in a single parallel bond ar- 
ray. For the Hertz distribution, the probability that a 
given bond has a radius greater than or equal to r, B > (r), 
is 



B>(r) 



(6) 



Then the radius distribution of the smallest bond from 
among w, S w (r), is given by ](| 



S w (r) = wb(r)B > (r) w - 1 =2awr e - m 



(7) 



The first equality expresses the fact that one of the w 
bonds has (smallest) radius r, with probability b(r), and 
the other w — 1 must have radii larger than r. 
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From the basic connections between permeability, pore 
radius, and time scale (Eqs. dJ) and (||)), and the fact 
that the clogging time is dominated by t\, we deduce 



k(w) 
wk(1) 



(8) 



while the clogging time distribution, P W (T), is directly 
related to the smallest bond radius distribution through 
P W {T) dT = S w (r) dr. Using Eqs. (g) and (|), we obtain 
the basic result (independent of system length) 



Pw(T) 



T 3/2 



o/T 1 ' 2 



(9) 



The power law should apply in the time range w 2 < t < 
w 2 N 2 . The former corresponds to the size of the typical 
smallest pore in a single realization, s/y/w, while the lat- 
ter corresponds to the smallest pore from among TV real- 
izations of the system, s/yNw. The essentially singular 
short-time cutoff arises from those realizations where the 
smallest bond happens to be anomalously large. Other 
coincident particle and bond radius distributions lead to 
qualitatively similar forms for P W (T), but with a differ- 
ent exponent value. For example, if S w (r) ~ as r — > 0, 
then the long-time exponent in P W (T) is — (jjl + 5)/4. 

We test our predictions by Monte Carlo simulations 
of the motion and trapping of suspended particles in a 
lattice porous medium. Each bond (with unit length) 
corresponds to a pore and the sites represent pore junc- 
tions. The algorithm is event-driven to bypass the phys- 
ical slowing down of the flow as clogging is approached. 
At an intermediate stage, there are a finite number of 
particles in the network, consistent with a unit pressure 
gradient and a constant-density suspension. Particles are 
defined to be always at lattice sites and each (including 
newly injected particles) evolves by either: (a) moving 
to the next "downstream" site, (b) blocking the down- 
stream bond it has entered, or (c) remaining stationary. 
The probability of each of these possibilities is defined 
so that the process corresponds to particles moving, on 
average, at the local velocity. 

To achieve this, each particle at a junction selects an 
outgoing bond i, with entrance probability equal to the 
fractional flux into this bond. In an update, a particle at- 
tempts to move to the next junction, via the pre-selected 
bond, with probability proportional to Vi/v max , where Vi 
is the local bond velocity and v max is the instantaneous 
largest particle velocity in the system. If the attempt 
occurs, the particle either moves to the next downstream 
site, or, if the particle is too large, blocks the bond. After 
a single update of all particles, the time is incremented 
by At — l/u m ax to ensure that, on average, a particle 
moving along bond i has speed Uj. 

After all particles undergo their move attempts, Af oc 
cj> At new particles are injected at the upstream end of 
the network, where <fi is the overall fluid flux. This en- 
sures the correct absolute velocity for each particle. Since 
the time increment systematically increases as bonds get 



blocked, the particle injection rate progressively slows. 
Our approach is distinct from earlier multiparticle simu- 
lations of filtration, where the injection rate was decou- 
pled from the overall flow p| . 
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FIG. 2. (a) Distribution of clogging times on the bubble 
model (continuous curve) and on a square lattice (data points) 
for system width 20 and length 10. Time units correspond to 
an initial injection rate of 0.05 particles per unit width per 
unit time. The straight line has slope —3/2. (b) The scaled 
clogging time distribution for square lattices of length 10 and 
widths 20, 40, 60, 100, and 200 (right to left). The data have 
been smoothed over a 12-point range. The heavy dashed line 
is the prediction of Eq. (g) . 

In Fig. ^(a), we plot the clogging time distribution 
from simulations on the square lattice and the bubble 
model of the same spatial extent. The latter is a series of 
parallel bond arrays with perfect mixing at each junction; 
this idealized system accounts for geometrical aspects of 
clogging [Q. The agreement between the two distribu- 
tions is remarkably good, suggesting that the schematic 
evolution proposed in Fig. [I] is quantitatively correct. 
The tails of the two distributions in Fig. ^(a) are well-fit 
by a power law with exponent —3/2. As a further test 
of Eq. ([)]), we plot, in scaled units, data for the clog- 
ging time distribution for systems of various widths w 
(Fig. ||(b)) and compare to the functional form in Eq. (||). 
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An important consequence of the power-law tail in the 
clogging time distribution is a transition in the corre- 
sponding moments. These are 

/>oo t>N 2 w 2 

(t k ) = tP w (t)dtn wt k ~ 3 ^ 2 dt, (10) 

JO Jw 2 

where the approximation replaces the short-time cutoff 
at t ps w 2 in Eq. (||) by the lower limit in the integral. 
Consequently 

( w 2 N 2-l/k fc>1 / 2; 

M k (w) = (i fc ) 1/fc ~ { w 2 {lnN) 2 k = 1/2; (11) 
[w 2 fc<l/2. 

Thus the mean clogging time diverges when the number 
of realizations is infinite, leading to large sample to sam- 
ple fluctuations in the clogging time, while the moments 
Mfe with k < 1/2 should be well behaved. The transi- 
tions in the moments is illustrated in Fig. (j^), where 
behaves erratically for k — 1 and 1/2, but then appears 
to grow as w 2 for k < 1/3. 
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FIG. 3. Width dependence of the moments (t k ) 1/k for 
k = 1,1/2,1/3,1/4, and 1/6 (top to bottom). The dashed 
line corresponds to a w 2 dependence. 

The behavior of Mk raises the question of when is a 
filter no longer useful. Waiting until complete clogging is 
impractical because of large fluctuations in the clogging 
time and eventual poor filter performance. It would be 
more useful to operate a filter only until the permeability 
decays to a (small) fraction of its initial value, such that 
reasonable flow and trapping efficiency are maintained, 
while minimizing fluctuations in this threshold time. We 
provide some preliminary results about these questions. 
Simulations indicate that the time to reach permeabil- 
ity fraction /, t(f), is proportional to u; 2 / 3 / -1 / 2 for 
10~ 10 < / < 10~ 3 , with the distribution of these frac- 
tional times progressively broadening for decreasing / 
(Fig. ^). Thus waiting until the permeability decays to 
a fixed fraction can provide at least a reliable criterion 
for when a filter should be discarded. A consequence 
of t(f) cx w 2 / 3 / -1 / 2 is that the permeability decays as 



1/t 2 over a wide range. Because of this rapid decrease, 
a practical filter should also have pores typically larger 
than particles to have a reasonable lifetime. These issues 
and understanding the efficiency of a filter throughout its 
useful life are currently under investigation. 




FIG. 4. Probability that a filter with w — 200 reaches 
permeability fraction / at time tj, for / = 4~™, with 
n — 2,4, 6, 8, 12, 16, together with the time distribution until 
complete clogging (left to right). 
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